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ABSTRACT 

We present HST /STIS observations of the optical counterpart (OT) of the 
7-ray burster GRB 000301C obtained on 2000 March 6, five days after the burst. 
CCD clear aperture imaging reveals a i? ~ 21.50 ± 0.15 source with no apparent 
host galaxy. An 8000 s, 1150 < A/A < 3300 NUV-MAMA prism spectrum 
shows a flat or slightly rising continuum (in fx) between 2800 and 3300 A, with 
a mean flux 8.7ttl ± 2.6 10-18 

ergs s 1 cm ^ A , and a sharp break centered 
at 2797 ± 25 A. We interpret this as the H I Lyman break at ^ = 2.067 ± 0.025 
indicating the presence of a cloud with a H I column density log (A^hi cm^) > 18 
on the hne-of-sight to the OT. This measured redshift is conservatively a lower 
hmit to the GRB redshift. However, as all other GRBs which have deep HST 
images appear to lie on the stellar field of a host galaxy, and as the large H I 
column density measured here and in later ground-based observations is unlikely 
on a random linc-of-sight, we believe we arc probably seeing absorption from 
H I in the host galaxy. In any case, this represents the largest direct redshift 
determination of a 7-ray burster to date. Our data are compatible with an OT 
spectrum represented by a power-law with an intrinsic index a — 1.2 {f^, oc i/"") 
and no extinction in the host galaxy or with a — 0.5 and extinction by SMC-like 
dust in the OT rest-frame with Ay — 0.15. The large Ahi and the lack of a 
detected host are similar to the situation for damped Ly-a absorbers at 2; > 2. 

Subject headings: gamma rays: bursts 

1. Introduction 

Gamma-ray bursts (hereafter, GRB) remain mysterious nearly 30 years after their first 
detection (Klebesadel, Strong, & Olson 1973), although some progress has been made in their 
understanding. The large sample of events collected by the BATSE experiment on-board 
the Compton Gamma-Ray Observatory shows an isotropic angular distribution (Meegan et 
al. 1992). However, their intensity distribution shows fewer weak bursts than expected from 
a homogeneous distributions of sources in a Euclidean space. It was therefore thought that 
GRBs were of extragalactic origin (Meegan et al. 1992). 



* Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract No. NAS5-26555. 
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Rapidly distributed arcminute localizations derived from X-ray instruments such as 
those on BeppoSAX (Costa et al. 1997) and, more recently, RXTE (Smith et al. 1999) 
satellites have revealed the existence of fading afterglows and enabled efficient searches for 
GRB counterparts to be carried out at other wavelengths. The optical transient (hereafter, 
OT) to GRB 970228 was the first to be discovered (van Paradijs et al. 1997). Its location 
angularly close to a faint diffuse object, possibly a galaxy, suggested that they are associated, 
providing further evidence that GRBs he at cosmological distances (van Paradijs et al. 1997; 
Sahu et al. 1997, later, Pruchter et al. (1999a) showed that the diffuse object is indeed a 
V — 25.8 galaxy). This hypothesis was soon confirmed when a direct spectrum of the OT 
to GRB 970508 revealed absorption fines at the same redshift as the underlying z — 0.835 
galaxy (Metzger et al. 1997; Bloom et al. 1998). Since then, OTs have been found for roughly 
half of the GRBs for which an X-ray afterglow has been detected and most of them appear 
to be associated with faint, z > 0.4 galaxies (Hogg & Pruchter 1999; Kulkarni et al. 2000). 
In particular, the probable host galaxy of GRB 971214 has a redshift of ^ = 3.42 (Kulkarni 
et al. 1998). 

Most popular theories describing the GRB phenomenon suggest that GRB activity will 
be closely tied in time to episodes of massive star formation. The collapse of massive stars 
(Paczyhski 1998) requires an active or recently active star forming host, while most neutron 
star - neutron star or neutron star - black hole binary mergers (Paczyhski 1991; Narayan, 
Paczynski, & Piran 1992) occur not long after star formation. The evidence for such an 
environment is already suggestive, as a number of host galaxies show strong emission lines 
associated with star formation (Metzger et al. 1997; Djorgovski et al. 1998; Kulkarni et al. 
1998; Bloom et al. 1999; Vreeswijk et al. 2000), while direct imaging with HST/STIS and 
NICMOS indicates that they are unusually blue with V — H < 1 whereas most galaxies in 
the Hubble Deep Field have V — H > 1 (Pruchter et al. 1999b). We thus expect that the 
spectra of OTs will show evidence that the events take place in environments of massive star 
formation: large hydrogen column density, large extinction by dust or even absorption by 
molecular hydrogen (Draine 2000). All lead to features in the rest-frame UV region. 

We have thus embarked on a program to obtain, as targets of opportunity, HST /STIS 
near ultraviolet spectra of GRB OTs. Although Bloom et al. (1997) have suggested using 
grating modes of STIS for studying GRBs, we chose to use the prism as it is the most efficient 
available spectral element due to its low dispersion and broad spectral coverage. Combined 
with the near ultraviolet, photon-counting multianode microchannel array (hereafter, NUV- 
MAMA) detector, the prism mode provides spectral coverage over the wavelength range 
1150 < A/A < 3300. However, its highly non-linear dispersion makes the wavelength cali- 
bration at the red end (2800 A/A < 3300) especially critical. 
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Here we present HST /STIS clear aperture CCD imaging and the first ultraviolet spec- 
trum of a GRB optical counterpart, GRB 000301C. This burst was detected at O'^Sl^ST* 
(UTC) on 2000 March 1 by the ASM instrument aboard the RXTE sateUite, as well as by 
Ulysses and NEAR (Smith et al. 2000). Its 7-ray light curve showed a common, simple 
shape: singly peaked, fast-rise followed by a slow decay, lasting a total of 10 seconds, which 
qualifies it as a short/intermediate duration burst (Jensen et al. 2000). The resulting com- 
posite localization of area ~ 50 arcmin^ was imaged on 2000 March 3.14-3.28 (UTC), by 
Fynbo et al. (2000a) with the Nordic Optical Telescope (hereafter, NOT) which revealed a 
then R — 20.09 ± 0.04, blue optical counterpart (Jensen et al. 2000). 

All exposures were taken on 2000 March 6 (UTC) during a single, 5 orbit visit. The 
first orbit was dedicated to imaging, acquisition and cahbration. Four 2000 s NUV-MAMA 
prism spectra were obtained during the subsequent four orbits. 

2. Imaging 

The field of GRB 000301C was observed using the STIS 50CCD, clear aperture mode. 
In this setting no filter is interposed between the CCD and the sky, and the bandpass is 
determined entirely by the reflectivity of the optics and the response of the CCD. Three 
exposures of 480 s each were taken in a diagonal dither pattern and combined using the 
Drizzle algorithm and associated techniques (Pruchter & Hook 1999; Pruchter et al. 1997). 
The drizzled output image, presented in Figure 1, has pixels of one-half the linear dimensions 
of the original images, or 0'.'025355 ± 0.000035, since the pixel scale of the CCD is 0'.'05071 ± 
0.00007 (Malumuth & Bowers 1997). 

At the position of the optical transient of GRB 000301C we find an unresolved object 
(FWHM 0'.'087) whose appearance is consistent with the STIS CCD Point Spread Function 
(hereafter, PSF). The mean count-rate of detected photons within 20 CCD pixels (containing 
more than 99.6% of the fiux (Leitherer et al. 2000)) from the OT image centroid is 56.0 counts 
s~^ on mean 2000 March 6.22 (UTC). Assuming a power law spectrum with a = —0.8 and a 
foreground Galactic extinction with Ay = 0.16, as determined from the dust extinction map 
given by Schlegel, Finkbeiner & Davis (1998), we obtain a magnitude oi R — 21.5 ± 0.15, 
where the uncertainty is dominated by the fact that the STIS CCD clear aperture mode is 
not completely calibrated. 

Any underlying galaxy would have to have i? > 24, or be unusually compact for it to 
avoid detection in our data. We therefore beheve that the apparent fiattening - the decrease 
of the OT luminosity fading rate - of the light curve in the R-hand reported by many 



-6- 



observers before this image was taken (cf. Garnavich et al. (2000); Halpern et al. (2000a); 
Bernabei et al. (2000); Rhoads & Fruchter (2001) and references therein) is not due to an 
underlying host galaxy. Instead, it may either be due to the jet itself (Bergcr ct al. 2000) or 
to micro-lensing (Garnavich, Loeb, & Stanek 2000). The diffuse emission to the northwest 
of the OT is the galaxy previously identified as a possible host (Rhoads & Fruchter 2000). 
From our STIS image, we derived an AB magnitude 25.1 ± 0.2 in the STIS clear aperture 
passband. The lack of emission between the two and the good fit between the OT and the 
STIS PSF suggest that the OT host galaxy has not yet been detected. 



3. Spectroscopy 

3.1. Acquisition and reduction 

In this section, we describe the acquisition and the reduction procedures which required 
special attention. The reduction was performed using IDL and the GSFG version of GAL- 
STIS (Lindler 1999). 

Each TIME-TAG, 2000 s spectrum was obtained with the prism using the 2125 A 
setting - so that the spectral region corresponding to wavelengths ~ 2125 A falls onto the 
central region of the NUV-MAMA detector - and the 52" x 0.5" slit. After each scientific 
spectrum an automatic wavelength calibration exposure was taken through the narrower 
52" X 0.05" slit. 

We used the it! = 18.05 star located ~ 5.7" west to the OT (cf. Fig. 1) as an offset star. 
The offset values were determined from the 3 x 900 s i?-band NOT image. Acquisition took 
place after the STIS CCD imagery and just before the NUV-MAMA prism spectroscopy. 
This order allowed us to obtain the CCD images independently of the success of the spectral 
acquisition. However, as the acquisition produces a small change in the location of the 
object on the detector, an additional offset has to be taken into account for the zero-point 
determination of the wavelength cahbration (see below) . 

Each scientific raw dataset displays a spectrum whose spatial profile can be modeled by 
a Gaussian with a mean FWHM of 2.5 pixels, or 0.072", as expected from a point source. 
There is no significant variation in the count rate between the different datasets, as well as 
no significant variation in the count-rate during the course of an orbit. 

The sum of the sky and dark backgrounds at the location of the OT spectrum was 
determined for each individual exposure. It was obtained by interpolation based on a fifth 
degree polynomial fitted independently to each column of the image, avoiding a 20 pixel 
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wide region centered on the object spectrum. Extraction was performed by summing up 7 
rows centered on the peak of the 2-D spectrum. 

Each extracted spectrum was then wavelength cahbrated. This step is especially critical, 
so we describe it here in more detail. 

The prism dispersion relation A = A(X) gives the wavelength A corresponding to the 
pixel^ number X in the extracted spectrum. It is well modelled by a relation 

Due to the strong blending of the hues especially at the red end of the spectrum, no deter- 
mination of the dispersion solution (i.e., the values of each coefficient Oj) is feasible on-orbit. 
Instead, the values of Oj were obtained during the pre-launch, thermal vacuum testing of 
STIS, using an external platinum discharge lamp fed through a vacuum monochromator, the 
52" X 0.05" slit and the same 2125 A setting as for the observations. The residuals to the 
fit were then measured to have an rms equal to 0.16 pixels. 

The quantity Xq is, in pixels, the separation along the dispersion direction of the pro- 
jected location of the object on the detector from the projected location of the center of 
the 52" X 0.05" slit used for the determination of the dispersion solution. Xq determines the 
zero-point of the wavelength calibration for each exposure and can be represented as the sum 
of three quantities: (1) A^row, the difference in pixels between the angular separations along 
the dispersion direction between the OT and the offset star as determined on the STIS CCD 
and NOT images; (2) Ax, the difference in pixels between the projection on the detector of 
the center of the slit used for the wavelength cahbration relative to its location during the 
pre-launch calibration; (3) A/9, the difference in pixels between the location of the center 
of the slit used during the scientific exposures (52" x 0.5") and that used for the automatic 
wavelength calibration after each scientific exposure (52" x 0.05"). In the following, we 
describe how each of these three quantities can be determined. 

The acquisition process usually centers the object in the slit used for the scientific 
exposure. However, the nearby R = 18.05 bright star was considered too close to the OT 
and could have interfered with the normal procedure. Instead, it was used as an offset star. 
Consequently, the acquisition process (a) changed the spacecraft orientation so that the 
offset star is centered on the slit by an on-board centroid algorithm, and then (b) modified 
the spacecraft pointing to take into account the offsets in right ascension and declination 



^Throughout this paper, a pixel refers to a native, 1024x 1024 format MAMA pixel, not a 'high-resolution' 
pixel (cf. Woodgate et al. (1998); Kimble et al. (1998)). 
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between the OT and the offset star as determined on the NOT image. However, the STIS 
image provides more accurate, but slightly different values for these offsets, so that we can 
determine exactly how far the object actually is from the slit center. A^j-ow represents the 
difference along the dispersion direction between the two offset determinations. Note that 
we also checked that the on-board algorithm provided a centroid location for the offset star 
on the acquisition image consistent, within 0.1 CCD pixel (0.2 MAMA pixel), with the one 
used to determine the offset between the star and the OT on the STIS CCD image. 

Ideally, Ax is best determined for each scientific dataset by cross-correlating its asso- 
ciated wavelength calibration exposure with a prc-launch wavelength calibration exposure. 
Unfortunately, the internal lamps were not operated during the prc-launch tests at the same 
current setting as during the on-board calibration, so that the shape of the lamp spectra 
were too different to be useful. In order to mitigate this effect, we built a template of the 
prism wavelength calibration lamp at the same setting as used during our observations based 
on G140L, G230L and G430L wavelength calibration exposures and the relevant sensitiv- 
ity curves. However, even this template does not show exactly the same shape at the red 
end as the wavelength calibration spectrum. The value of Ax was finally obtained by cross- 
correlating a 61-pixel long region (along the dispersion direction) of each spectrum where the 
dispersion is relatively large. As a check, we measured the center of the Lyman-a geo-coronal 
line determined by the middle position between its edges, and found that it falls within 1/6 
of a pixel of the expected location. Finally, shifts due to thermal or flexure causes were 
measured to be between 0.02 and 0.2 pixels from orbit to orbit, and thus probably less than 
0.1 pixel during an exposure. 

Finally, regular on-board STIS calibrations enable us to determine A/3 to better than 
1/20 CCD pixel, (1/10 MAMA pixel). 

Consequently, a conservative (2(t) estimate of the error on the zero-point of the wave- 
length calibration is 0.5 MAMA pixels, or 0.23, 1.4, 5.2, 12, 22 and 29 A at A = 1200, 1500, 
2000, 2500, 3000 and 3300 A, respectively. 

Since the dispersion is small at ~ 2900 A (the location of the peak), the effect of the 
telescope and instrument line spread function (hereafter, LSF) must be taken into account 
to obtain a correct estimate of the count rate per pixel. We, therefore, deconvolved the 
observed spectrum using the Lucy- Richardson method. Unfortunately, the blending of lines 
in the wavelength calibration spectrum prevents its use to create a LSF. In the absence of 
suitable calibration frames taking into account the effect of the Point Spread Function of 
the telescope itself, we use a LSF based on the shape perpendicular to the dispersion of the 
2-D spectrum of the flux standard star HS 2027+0651 at about 2800 A. Ray-tracing studies 
indicate that this is a good approximation of the LSF along the dispersion direction (C.W. 
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Bowers, private communication). This star was observed in the prism mode during a STIS 
cahbration program. The immediate consequence is that the corrected (de-convolved) mean 
level of the count rate is larger by about 25% between 2800 and 3300 A compared to the 
un-corrected one. 

Before converting the count rate to flux, we have checked the wavelength calibration of 
the prism spectrum of HS 2027+0651 and deconvolved it before building a new sensitivity 
curve. Once converted to flux, the 4 individual spectra were rebinned to a same set of 
wavelengths using bi-linear interpolation and averaged. The final, flux calibrated spectrum 
of GRB 000301C is presented in Figure 2a. 

3.2. Simulations 

Since Lucy- Richardson deconvolution of a low S/N spectrum actually decreases its S/N, 
we prefer to compare the spectrum in total observed counts per pixel in wavelength with 
simulated spectra. These spectra were built using an input spectrum expressed in units of 
flux, multiplied by the sensitivity curve, integrated over the spectral range of each pixel, 
and finally convolved with the LSF. The output spectrum expressed in count-rate is then 
multiplied by the exposure time. 

In this paper, we only consider input spectra behaving as power-law of the form /j, oc 
as expected for GRB afterglows (Meszaros, Rees & Papathanassiou 1994; Sari, Piran & 
Narayan 1998), modified to take into account the Galactic extinction towards GRB 000301C, 
the Lyman-a forest absorption lines, possible absorption due to H I, H2 as well as extinction 
by dust located in the host galaxy. 

For the Galactic extinction in the direction towards GRB 00030 IC, we used the value 
Ay — 0.16 as determined from the dust extinction map given by Schlegel, Finkbeiner & 
Davis (1998) and the analytical expression of the Galactic extinction curve given by Pei 
(1992). In order to model the Lyman-a forest absorption spectrum, we first noted that 
the likely redshift of GRB 000301C (cf. below) falls in a range (1.5 < z < 2.3) where the 
Lyman-a forest is still poorly studied. So far, the Hubble Deep Field South (HDF) quasar is 
the only one in that range whose absorption spectrum has been observed at resolution high 
enough to determine the H I column density and Doppler parameter by Voigt profile fitting. 
We thus used the H I absorbers listed by Savaglio et al. (1999) to estimate an absorption 
spectrum due to their respective Lyman series (Ly-o; to Ly-12) lines. However, we discarded 
the fines associated with the z — 1.942 Lyman-fimit system seen in the HDF quasar as it is 
a single event whose occurrence is dominated by small number statistics. Naturally, we also 
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eliminated the lines due to clouds whose redshift is larger than the assumed GRB redshift 
in the simulation. 

Absorption due to H I located in the host galaxy takes into account Lyman series lines, 
modeled by Voigt profiles, as well as continuous absorption. The redshift of the host galaxy 
the H I column density A^hi and the Doppler parameter h can be easily adjusted. The 
modelling of possible absorption due to H2 located in the host galaxy follows Draine (2000). 

Finally, the possible extinction due to dust in the host galaxy, characterized by the 
rest-frame Ay extinction and the type of dust (Small Magellanic Cloud, Large Magellanic 
Cloud or Milky Way) also uses the expression given by Pei (1992). 

3.3. Results 

Figure 2a reveals a flat or slightly rising continuum spectrum (in fx) over the range from 
2800 A to the sensitivity limited red end of the spectrum (3300 A). No useful constraint 
on the spectral slope can be derived from the spectrum alone. The mean flux is measured 
to be 8.7+?:^ ± 2.6 lO'^^ ergs s ^ cm ^ A , where the first error values come from the 
uncertainties in the wavelength calibration combined with the strong slope of the sensitivity 
curve at ~ 2900 A; the second error value is the standard deviation of the flux calibrated 
error array over the relevant pixels. However, a sharp break can be seen centered at 2797±25 
A(2(t). There is no significant fiux recovery at the blue end of the spectrum. 

3.4. Redshift of GRB 000301C 

The observed break is best interpreted as the onset of continuous absorption below the 
H I Lyman break due to the presence of a cloud with a large neutral hydrogen column 
density, located at z = 2.067 ± 0.025 (2(t) in the line-of-sight to the OT. From the fact that 
the observed fiux is zero within the error bars blueward of the break, we determine that 
log(AHi cm^) > 18.0. The decreased sensitivity at the blue end combined with the possible 
presence of other relatively high column density clouds in the Lyman-cc forest does not allow 
us to set a higher lower hmit. An upper limit log (Ahi cm^) < 23.3 can be set from the 
absence of a strong feature at the expected location of the Ly-/3 line at ~ 3100 A. 

Our redshift determination is confirmed by the presence of weakly or marginally detected 
lines, compatible with Fe II, Mg II and other low- ionization species at a redshift z — 2.0335 ± 
0.0003 in a Keck spectrum obtained by Castro et al. (2000). More importantly, an ESO/VLT 
spectrum covering the wavelength range from 3600 A to 8220 A, was obtained by Jensen et 
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al. (2000). It reveals a damped Ly-a with a large, although uncertain log A^hi = 21.2 ± 0.5, 
within the range allowed by our spectrum, as well as some associated C IV, Fe II and possibly 
other low-ionization lines. These lines appear ai z = 2.0404 ± 0.0008, a significantly larger 
value than Castro et al. (2000) and in reasonable agreement with our value. 

The value z — 2.067 ± 0.025 corresponds to the largest direct redshift measurement 
of the OT to a GRB so far. If we assume that the H I cloud is associated with the OT 
host galaxy, then the above value is the redshift of the OT. We note however that it is 
conservatively a lower limit to the OT redshift. Indeed, if the object giving rise to the GRB 
is located outside of a galactic halo or, more generally, any H I cloud, all the absorption lines 
detected in the OT spectrum would be due to intervening systems. 

Can we set an upper limit to the GRB redshift? Although the ESO/VLT spectrum 
extends down to 3600 A, its low signal-to-noise ratio at the blue end cannot allow us to 
exclude the presence of low column density Ly-a forest lines between 3600 A and ~ 4000 
A. This wavelength corresponds to a firm upper limit oi z — 2.3 to the GRB redshift. This 
value also corresponds to the maximal redshift for which our STIS spectrum could reveal 
Lyman-limit systems which have log (A^hi cm^) > 17. Using the number density of Lyman- 
limit systems seen in QSO spectra (Stengler-Larrea et al. 1995) in the relevant redshift 
range, we estimate that the probability of observing no additional Lyman-limit system over 
the range 2.067 < z < 2.3 to be 0.68. This value is too small to exclude a larger redshift 
than z = 2.067 ± 0.025. Similarly, the probability that a random linc-of-sight covering the 
redshift range < 2; < 2.3 would encounter at least one damped Ly-a system is ~ 0.25 
based on the number density of damped Ly-a sytem determined by Rao & Turnshek (2000). 
This value is too large to ascertain that the redshift of the high-column density system 
corresponds to the redshift of GRB 000301C. In other words, our HST/STIS spectrum and 
the ESO/VLT one only allow us to hmit with certainty the redshift range of the OT to 
be 2.067 < z < 2.3. However, we can argue in favor of the z — 2.067 ± 0.025 value for 
the redshift of GRB 000301C. Indeed, on one hand, damped Ly-a systems are thought 
to be mainly caused by merging protogalactic clumps hosted by collapsed DM halos and, 
consequently, forming stars as the progenitors of normal present-day galaxies (Gardner et al. 
1997; Haehnelt, Steinmetz & Ranch 1998, 2000). On the other hand, most OTs discovered 
so far are found to be associated with galaxies (Pruchter et al. 1999b), many of which show 
signs of significant star formation (Metzger et al. 1997; Djorgovski et al. 1998; Kulkarni et al. 
1998; Pruchter et al. 1999a,b; Bloom et al. 1999; Vreeswijk et al. 2000) for which, presumably, 
a large reservoir of gas is available. In summary, the unusually large H I column density and 
the fact that most OT host galaxies found so far appear to actively form stars lead us to 
believe that the absorption seen in the spectrum of GRB 000301C is most likely due to the 
host galaxy. 



3.5. Low extinction towards GRB 000301C 



When combined with the contemporaneous STIS CCD imaging, the STIS spectrum 
also allows us to constrain the extinction taking place in the z = 2.067 H I cloud towards 
GRB 000301C. Figure 2b compares the observed total number of counts to a simulation 
characterized by an intrinsic power-law spectrum with index a — 0.5, and some extinction 
in the host galaxy by SMC-like dust {Ay — 0.15). However, in the wavelength range of the 
STIS spectrum, such a model is indistinguishable from one defined by an intrinsic power-law 
spectrum with index a = 1.2 (/,^ oc and no extinction in the host galaxy. These two 
models not only are able to reproduce the appearance of our spectrum without the need of 
any molecular hydrogen, but also conform to the slope of the optical spectrum (Halpern et 
al. 2000b; Feng, Wang & Wheeler 2000; Jensen et al. 2000) within its uncertainty. They 
also match the count-rate of the STIS CCD image, within the 15% accuracy of the IDL 
routine SIM_STIS (Plait, private communication), which predicts the expected count-rate 
for the different STIS modes based on the individual optical components. However, if a 
steep spectrum with a = 1.2 is compatible with our STIS imaging and spectroscopic data 
alone, a more shallow one {a ~ 1) is favored when ground-based observations with lower 
uncertainties in optical band magnitudes than the STIS-derived R magnitude are taken into 
account (cf. Halpern et al. (2000b); Rhoads & Fruchter (2001)). Similarly, a < 1, models 
with low extinction Ay < 0.1 {Ay < 0.05) due to LMC- (Milky Way-) like dust are also 
consistent with the STIS data alone. They predict that the 2175 A feature (cf. e.g. Blanco, 
Fonti & Orofino (1996)) should be detectable and shifted to ~ 6570 A. However, its presence 
seems to be ruled out by existing optical spectra and photometry (Rhoads & Fruchter 2001; 
Jensen et al. 2000). Therefore a model with a relatively small amount (Ay < 0.1) of SMC-like 
dust appears to best fit the full UV/optical data. 



4. Discussion and conclusion 

The HST/STIS CCD observations of GRB 000301C revealed aR = 21.50 ±0.15 source 
on 2000 March 6.2 (UTC) with no trace of a host galaxy. The NUV-MAMA prism spectrum 
presents a relatively flat spectrum (in fx) between 2800 and 3300 A, with a mean flux 
8.7+?;g±2.6 10"^^ ergs s~^ cm^^ A~^ and a sharp break centered at 2797 ±25 A, interpreted 
as a H I Lyman break. It indicates that the OT arises in or beyond an absorber ai z = 
2.067 ± 0.025 similar to the ones causing the high column density systems detected in QSO 
spectra. The redshift and the high H I column density are confirmed by the observed 
logA^Hi = 21.2 ± 0.5 damped Lj-a, Fe II, C IV and other lines observed in the ESO/VLT 
spectrum (Jensen et al. 2000), as well as the Mg II and Fe II lines in the Keck spectrum 
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(Castro et al. 2000). The lack of a detected host galaxy and the large H I column density 
system arc similar to the situation for damped Ly-a systems at z > 2 (e.g. Lowcnthal et al. 
1995; Fynbo, Moller & Warren 1999), which are thought to be mainly caused by merging 
protogalactic clumps hosted by collapsed DM halos and consequently probe the progenitors 
of normal present-day galaxies (Gardner et al. 1997; Haehnelt, Steinmetz & Ranch 1998, 
2000). Combined with the fact that other OTs are often found in star forming galaxies, 
presumably with large H I content, this strongly suggests that z — 2.067 ± 0.025 is actually 
the redshift of GRB 000301C. If this is the case, our spectrum as well as the ESO/VLT and 
the Keck ones suggest that the hne-of-sight to the OT crosses the interstellar medium of its 
host galaxy. 

Our data are compatible with no or little extinction. The absolute magnitude of the 
OT at the time of the discovery image, 41.5 hours after the burst, was therefore Mr = —26.1 
for ifo = 65 km s~^ Mpc~^, Q = 0.3 and A = 0.7. The low extinction would argue against a 
naive interpretation that 7-ray bursts originate in a star forming region. Instead, one might 
imagine this favors the hypothesis that GRB 000301C results from a neutron star - neutron 
star or neutron star - black hole binary merger (Paczyhski 1991; Narayan, Paczyhski, & Piran 
1992). In this case, this relatively short duration event could have taken place at a significant 
distance from the host galaxy some time after a significant burst of star formation. Narayan, 
Paczyhski, & Piran (1992) predict that the median distance between a binary merger and its 
host galaxy is ~ 50 kpc at the time of the burst. Consequently, it is not impossible that the 
nearby galaxy ~ 2" to the NE (Rhoads & Pruchter 2000; Kobayashi et al. 2000) is actually 
the host, although it is very red where as previous hosts tend to be blue in optical - infra-red 
colors: we derive a R — K' color of 4.9 ± 0.5 based on magnitudes measured by Rhoads & 
Pruchter (2000) and Veillet (2000). In particular, we cannot exclude that this galaxy lies 
at 2; ~ 2.067 on the basis of the existing photometry. If this is the case, its H I disk must 
extend over at least ~ 20kpc (Hq = 65 km s~^ Mpc^"*^, Qq = 0.3 and A = 0.7) to cover the 
line-of-sight of the OT. This value is similar to the separation between the line-of-sight to 
the QSO Q 2233+131 and the galaxy identified by Djorgovski et al. (1996) as responsible 
for the z — 3.150 damped Ly-a in its spectrum. Spectroscopy of the galaxy ~ 2" to the NE 
of GRB 000301C is necessary to settle the issue. 

However, the present data cannot exclude that the burst itself destroys dust as high- 
energy photons find their way out of a dense cloud (Waxman & Draine 2000). On the other 
hand, the time elapsed between the formation of the GRB progenitor and the GRB event 
itself could be long enough for the progenitor to leave regions of large extinction as is the 
case for core-collapse supernovae (type Ila, Ib/c), a significant fraction of which are also 
little affected by extinction (e.g. Schmidt et al. 1994). Finally, the low extinction could be 
possibly explained if the progenitor of GRB 000301C was located in a low metallicity galaxy 
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which has only recently started to produce stars. The host of GRB 000301C might therefore 
resemble the galaxies giving rise to damped Ly-a absorbers. 

Note: After this paper was submitted, the direct determination of the redshift for 
the optical counterparts of two other GRBs have been reported with redshifts similar or 
larger than the one of GRB 000301C: GRB 000926 at ^ = 2.066 (Fynbo et al. 2000b) and 
GRB 000131 at ^ = 4.50 (Andersen et al. 2000). 
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Fig. 1. — A 10"xlO" portion of the HST /STIS clear aperture CCD image of the field of 
GRB 00030 IC showing the OT on the left and the offset star used for acquisition on the 
right. North is up and East is to the left. The ring visible on the left of the offset star is a 
ghost due to internal reflections in the STIS CCD window. 
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Fig. 2. — (a) Top: Deconvolved, flux calibrated, UV spectrum of GRB 000301C. A break 
is clearly seen at A2797 A. If caused by the onset of Lyman continuous absorption due to 
H I gas associated with the host galaxy, the OT redshift is z = 2.067 ± 0.025. (b) Bottom: 
Comparison of the observed UV spectrum, expressed in observed total number of counts with 
simulated spectra. The dotted curve corresponds to a model with a — 0.5 and extinction by 
a SMC-hke dust with Ay — 0.15. A model with a — 1.2 and no extinction would present a 
spectrum indistinguishable in this wavelength range from the one shown here. 



